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Abstract 

We consider the Schrodinger operator on zigzag graphs with a periodic potential. The 

P^ i spectrum of this operator consists of an absolutely continuous part (intervals sepa- 

^^ ' rated by gaps) plus an infinite number of eigenvalues with infinite multiplicity. We 

,J^ ' describe all compactly supported eigenfunctions with the same eigenvalue. We define 

C^ , a Lyapunov function, which is analytic on some Riemann surface. On each sheet, the 

Lyapunov function has the same properties as in the scalar case, but it has branch 

points, which we call resonances. We prove that all resonances are real. We deter- 
^sq . mine the asymptotics of the periodic and anti-periodic spectrum and of the resonances 

^ I at high energy. We show that there exist two types of gaps: i) stable gaps, where 

the endpoints are periodic and anti-periodic eigenvalues, ii) unstable (resonance) gaps, 
^^ . where the endpoints are resonances (i.e., real branch points of the Lyapunov function). 

■^ I We obtain the following results from the inverse spectral theory: 1) we describe all 

finite gap potentials, 2) the mapping: potential - all eigenvalues is a real analytic 

f^ ■ isomorphism for some class of potentials. 

fH ! We apply all these results to quasi-lD models of zigzag single-well carbon nanotubes. 

-(— > 



1 Introduction and main results 



> 

S^ . Consider the Schrodinger operator H = —d + q with a periodic potential q on the zigzag 

H ! graph r^. The zigzag graph T^ is related to a zigzag nanotube T^. We shall discuss this at 

the end of this section. In order to define T^ we introduce the so-called honeycomb lattice 

r = Ua;!^^ C M^, cu = (n,j, fc) e Z X JI X Z, JI = {0, 1, 2}, where the edges F^ are given by 

r^ = {a; = r^ + te^, t G [0, 1]}, e^ = e^- G M^ r(„,i,fc) = ^(n,o,fc) + ^o, 

r(nfi,k) = r{n-i,2,k) = n(eo + ei) + k{ei + 62), cq = (0, 1), ei = ^{^y^, 1), 62 = 2*^^' ~^)' 

see Fig. [T] and |21 The edges F^^ of length 1 form regular hexagons. Each edge F^ has 
orientation given by e^j = e,- with the starting point r^^ G M^. We have the coordinate 
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Figure 1: Zigzag graph F^ for (a) A^ = 3; (b) A^ = 1. The fundamental domain Fq is marked 
by a bold line. 



X = r^ + te^^ and the local coordinate t G [0, 1]. We define an oriented zigzag graph F^ by 
identifying F„ .,■ ^ with T^j^k+N for all (n, j,k) G Z x J x Z. Here and below for simplicity we 
write F„ j-fc = T(^n,j,k) and fn,j,k = f(n,j,k) etc. Thus F^ is a topological space in the quotient 
topology. As example, the graphs F^ and F^ are illustrated in Fig. 1. We represent F^ by 
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Z/{NZ), J = {0,1,2}. 



In our paper we assume that A^ = 2771 + 1 for some integer m, ^ 0. For each function y 
on F^ we define a function y^ = y\r^,uj G Z. We identify each function y^ on F^ with 
a function on [0, 1] by using the local coordinate t G [0, 1] and define the Hilbert space 
L^(F^) = (ScogzL'^O^uj)- Introduce the space C(F^) of continuous functions on F^ and 
the Sobolev space W'^(T^) = {y, ^u^ezu'^ e L^(F^); y satisfies the Kirchhoff Boundary 
Conditions: y G C(F^) a7id the following identities hold 



-y'nfiM + l/n,l,fc(0) - l/^,2,fe-l(l) = 0, Z/Ul,0,fc(0) - Z/^,,l,fe(l) + Z/^,,2,fe(0) = 0. (1.1) 

for all {n,k) G Z x Z^r.} Condition (|l.ip means that the sum of derivatives of y at each 
vertex of F^ equals and the orientation of edges gives the sign ±. Our operator H on F^ 
acts in the Hilbert space L^(F^). We define our operator H by {Hy)^ = —y'^ + qy^, where 
y = {yui)u>GZ e T){H) = W'^{r'^), q G ^^(O, l). in the case g = 0, we denote the operator H 
by if°. The operator if° is self-adjoint, see [Cal], and so is H, see Sect. 3. 

Recall the needed properties of the Hill operator Hy = —y" + q{x)y on the real line 
with a periodic potential q{x + 1) = q{x),x G M. The spectrum of H is purely absolutely 
continuous and consists of intervals a„ = [\^_i, \~],77 ^ 1. These intervals are separated 




v^S^^Nf*^ 



r-1,0,1 I '--1,0,2 l''-l,0,3 ; 



'■-2,0,1 



Figure 2: The unrolled honeycomb lattice. The zigzag graph is marked by bold, N=3. The 
fundamental domain is marked by a double bold line. Each edge F^ has orientation given 
by e^ = ej,j = 0, 1, 2 with the starting point r^ G M^. 



by the gaps 7„ = (A~, AJ) of length |7„| ^0. If a gap 7„ is degenerate, i.e. |7„| = 0, 
then the corresponding segments a„, a„+i merge. For the equation —y" + q{x)y = \y on the 
real line we define the fundamental solutions "dlx, A) and <f{x, A),x G M satisfying ■(9(0, A) = 
ip'{0,\) = l,'i9'(0. A) = (/9(0, A) = 0. The corresponding monodromy matrix Ai and the 
Lyapunov function A are given by 



M{X) 



^(1,A) v^(l,A) 
^'(1,A) ^'(1,A) 



A(A) 



V9'(1,A)+^(1,A) 



AgC. 



:i.2) 



The sequence \q < \^ ^ \^ < is the spectrum of the equation —y" + qy = \y with 

periodic boundary conditions of period 2, that is y{x + 2) = y{x),x G M. Here equality 
A~ = A^ means that A^ is an eigenvalue of multiphcity 2. Note that A(A^) = (— l)", n ^ 1. 
The lowest eigenvalue Xq is simple, A(A[|') = 1, and the corresponding eigenfunction has 
period 1. The eigenfunctions corresponding to A^ have period 1 if n is even, and they are 
anti-periodic, that is y{x + 1) = —y{x), a; G M, if n is odd. The derivative of the Lyapunov 
function has a zero A„ in each "closed gap" [A~, A^], that is A'(A„) = 0. Let /i„,n ^ 1, be 
the spectrum of the problem —y" + qy = \y,y{0) = y{l) = (the Dirichlet spectrum), and 
let z/„, n ^ 0, be the spectrum of the problem —y" + qy = \y, y'{0) = y'{l) = (the Neumann 
spectrum). It is well-known that fin,t^n € [-^nj^n] ^^^ ^o ^ ^o- Moreover, a potential q is 

even, i.e., q G -^Len(O) 1) =^ 9 ^ -^^(0) 1) ■ ?(1 - x) = q{x),x G [0, 1] \ iff |7„| = |;U„ - z/„| for 



all n ^ 1. Define the set an = {fin, n ^ 1} and note that a^ = {X E C : (/?(!, A) = 0}. 



For simplicity we shall denote Ta,i C T^ by Fq,, for a = {n,j) G ^i = Z x J. Thus 
F^ = Ua^Zi^a, see Fig [T] In Theorem 11.11 we will show that H is unitarily equivalent to 
©^iffc, where the operator H^ acts in the Hilbert space L^(F^) and is given by 

{HkfU = -f'L + g/a, ®a^zX e ^'(r^), (i.3) 

and the vector function / = {fa)a(^Zi satisfies the Kirchhoff conditions: 

/n,0(l) = /n,l(0) = s'^fnM, /n+l,o(0) = /„,l(l) = /„,2(0), S = e^t, (1.4) 

-/;o(l) + /;i(0) - s'jl,il) = 0, /;+,,o(0) - /;,(1) + /^^(O) = 0. (1.5) 

We reduce the spectral problem on the graph F^ to some matrix problem on M. In order 
to describe this we define the fundamental subgraph Fq of F^ by Fq = U'^^qToj, see Fig. [TJ 
On F^, the group Z acts via p o F;'^ ^ = F^_,_pj, {n,j,p) E Zi x Z. Thus Fq is a fundamental 
domain associated with this group action of Z. 

For the operator Hk we construct the fundamental solutions Qk{x, A) = {Qk,a{x, \))aeZi, 
and $fc(x, A) = {^k,a{x, X))aeZi, {x, A) G M x C, which satisfy 

-fa + g/a = ^/«, the Kirchhoff Boundary Conditions (0,(1131), (1.6) 

efc,;3(0,A) = <l>l,^(0,A) = l, e',,^(0,A) = <l>,,^(0,A)=0, /3=(0,0). (1.7) 

We introduce the monodromy matrix 

which is determined by 0^, ^k on the fundamental domain Fq. 

There are two methods to study periodic differential operators. The direct integral 
analysis usually used for partial differential operators [ReS] gives general information about 
the spectrum, but no detailed results. The method of ordinary differential operators, based 
on the Floquet matrix analysis, gives detailed results. Note that there are a lot of open 
problems [BBK] even for the Schrodinger operator with periodic 2x2 matrix potentials on 
the real line . 

We introduce the monodromy matrix, similar to the case of Schrodinger operators with 
periodic matrix potentials on the real line, see [YS]. After this, using the approach from 
[BK],[BBK],[CK] we introduce the Laypunov functions and study the properties of this 
functions, similar to the case of the Schrodinger operator with a periodic matrix potential 
on the real line. This is a crucial point of our analysis. Here we essentially use the results and 
techniques from the papers [BK], [BBK],[CK]. The recent papers [BBK],[CK] are devoted 
to Schrodinger operators with periodic matrix potentials (a standard case) on the real line. 
Remark that Carlson |Caj studied the monodromy operator to analyze the Schrodinger 
operator on a product of graphs. His results do not cover our case. 

We formulate our first preliminary results. 



Theorem 1.1. (i) The operator H is unitarily equivalent to (B'^Hk, where the operator H^ 
is given by il.!^) - n~^) . 

(ii) For any A G C \ (Td and k G Z^v there exist unique fundamental solutions 0fc, $fc of the 
system / li.6|) , (|i.7| ). Moreover, each of the functions Qk{x, A), $fc(x, X),x G F^ is meromorphic 
in X E C\a£) and each matrix Aik{X) satisfies 

In particular, 

where Ck = cos^, s^ = sin^s = e'^^ . Furthermore, IZM-k^^^^ = TkTZAiTZ^^ is an entire 
matrix-valued function. In particular, the function Dkij, A) = det(7Vlfc(A) — r/2), where I2 
is the 2x2 identity matrix, is entire with respect to X,t E C 

Remark that in contrast to the Schrodinger operator with periodic matrix potential on 
the real line (see |YS. or |CKj ). the monodromy matrix Ai^ has poles at the points A G o"/), 
which are eigenvalues of Hk, see Theorem II .21 Such a phenomenon has already been observed 
e.g. in ISH, PH, |K^ - 

Define the subspace Hk{X) = {ip E T){Hk) : Hki^ = Xip} for A G app{Hk),k G Z^. If 
dim?-^fc(Ao) = 00 for some Aq G app{Hk), then we say that {Aq} is a fiat band. In Theorem 
11.21 we describe all fiat bands and supports of eigenfunctions (see Fig. E)). 

Theorem 1.2. Let (A, A;) E ctj:) x Z^v- Then 
(i) Each eigenf unction from 7ik{X) vanishes at all vertices ofT^. 

(ii) Let (p = (p{-, A) I [0,1] and put c = (p'{l, A). Define the function ijj'^^'^ = {^a )aeZi by : 
if 1] = 1 — s^c^ 7^ 0, then 

i^n] = 0>/or all n^O, -1, j E Z3, and ip^^l = rjip, i/j^°l = op, ip^^l = c^ip, 

iJ% = 0, ^i°ii = -s'c^, ij% = -^, (1.11) 

if rj = 0, then 

<d = 0, V-S = <^' ^oi = c<^' ^ni = 0' «^^ ^ + O'i' e Z3. (1.12) 



Then each ^*^"^ = i^n-rri,j)(m,j)eZi G Ti-kiX), n E Z, and each f E Hk{X) has the fi 



orm 



/ = E/«^^"^ f-=V^/f:f^ V^l (/.)nez€f. (1.13) 

Moreover, the mapping f — > {/„}nez is a linear isomorphism between TikiX) and C."^ . 
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Figure 3: The supports of the eigenfunction tjj^^^: (a) r/ = 1 — s''ip'{l, A)^ = 0; (b) 1] ^ 0. 
Let Tk,± be the eigenvalues of A4k,k G Z^. Using (jl.lOjl we deduce that 

(1.14) 
If A; = 0, then we introduce the standard Lyapunov function Aq = ^'^^° = |(to,+ + 7^), 
which is entire. If fc 7^ 0, then we define the Lyapunov functions A^-i- = \{Tk,± + :;r^) (see 
[BBK],[CK]) and using ()1.14|) we get A_k,± = Ak,±- Below we prove the following identities 

Ak,± = ^k±y/pk, ^k = Ao + sl, p,, = -^{cl-Q, Cfc = COS—, Sfc = sin— . (1.15) 
If g = 0, then we denote the corresponding functions by A°,p°, ... In particular, we have 

Introduce the two sheeted Riemann surface TZk (of infinite genus) defined by ^/pk■ The 
functions Ak,±, /c G m = {1, 2, .., m}, are the branches of A^ = ^^ + y/pk (see FiglH) on the 
Riemann surface TZk- We describe spectral properties of Hk in terms of A^. 

Theorem 1.3. (i) The Lyapunov functions A^ -t, k Em satisfy / ti.i,5j) . 
(a) For any k G "Zn the following identities hold: 

a{Hk) = (rpp{Hk)Uaac{Hk), (Tpp{Hk) = (Td, ^^^(i^fe) = {A G M : Afc(A) G [-1, 1]}. 

(1.17) 
(Hi) If Afe(A) G (—1, 1) for some A G M, A; = 0, .., ?7i and A is not a branch point of Afc(A), 
then A'^(A) ^ 0. 

Remark. 1) If we know Aq, then we determine all Ak,Pk,k G Z^r by (jl.isp . 2) If we 
know pj for some j G m, then by p.l5j) . we determine all Ak, pk,k G Zjy. 

Definition 1. A zero of pk,k Efn = {l,2,..,m} is called a resonance of H . 

Let A^2n ^^^ ^02n+ij''^ ^ bc the zeros of det(7Wo — h) and det(A1o + h) respectively 
(counted with multiplicity). In Theorem 11.41 we will show that the periodic eigenvalues A^2n 
and the anti-periodic eigenvalues A^2n+i satisfy 

^oiK,n) = (~1)") -^0,0 < ^0,1 < ^0,1 < ^0,2 ^ -^0,2 < "^0,3 < "^0,3 < "^0,4 ^ ^0,4 < •■•• (1-18) 




Figure 4: Typical Lyapunov functions and the zigzag nanotube. The function A3 is shown 
forgGLe,en(0,l),iV = 6. 

For Hq we introduce the spectral bands cro,n = [Ajj'^.i, Xq^] ci^nd the gaps 7o,n = {Xq^, A,!^), 
n ^ 1. Using the asymptotics of the fundamental solutions i), ip, we obtain 
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as |A| -^ 00. The following theorem describes the basic properties of Hq. 



Theorem 1.4. The function Aq is entire and has the following properties: 

(i) The function Aq has only real simple zeros \o^n,n ^ 1, which are separated by the simple 

zeros TjQ^n of Aq: 770,1 < Ao,i < 770,2 < Ao,2 < r7o,3 < ••• and satisfy 

5 
i^n.y^n e 7o,2n, Ao(Ao,2n) ^ 1, Ao(Ao,2„-i) ^ "-, f OT any n^l. (1.20) 

(a) The periodic and anti-periodic eigenvalues A^„, n ^ satisfy il.l^) and have asymptotics 
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' I 0(1) 



'1.211 



n 



where cf) = arcsin | e [0, |] and q-n = Jq q{t)e^'^'^^^dt,qsn = Inig„. Moreover, (— l)"'Ao(A) ^ 1 
for all A G [An„, X^J and Ao,„ G [An„, X^ ]. 



(m) \ln\ = I/in -^n\ iff In = 70,2n- 

(iv) 7o,2n C 7„ for all n ^ 1. Moreover, for fixed n ^ 1 we have 7o,2n = iff In = ^- 
(v) Ao'(Ai+2n) = -I for alln^OzjJqe L^JO, 1). 

Remark that the last theorem gives the asymptotics of the gap length. 

Let A^2n ^^^ '^^2n+i^^ ^ be the zeros of det(A^fc — I2) and det(A^fc + 12). Below we 
will show that the periodic eigenvalues X^2n ^^^ ^^^ anti-periodic eigenvalues X^2n+i satisfy 
the equations Aa:(A^^) = (—1)". In Theorem 11.51 we show that A^^ satisfy the equation 



Ao(Aj2n) = COS ^^, Ao(Aj2„+i) = -1, ke m, (1.22) 



27ik 
and the labeling is given by: each A^„ is double and 



■^0,0 < ^1,0 < ^2,0 < •■■ < -^m.O < ^0,1 < ^0,1 < ^»n,2 < ^m-1,2 < •■■ < ^0,2 ^ "^0,2 < ■■•) (1-23) 

Aj2n-i = Aj2n-i, {k,n)emxK (1.24) 

The periodic eigenvalues A^2n' ^ G m = {1, .., m}, ra ^ (i.e., Afc(A^2n) — ^ ) satisfy 



\,2n = (^^ ± <^fc) +9oH , 0fc=2a^ccos ^[O'o] "'^ Tl ^ oo. (1.25) 

Let r^^, /c G m, n ^ be the zeros of p^. Below we will show that r^^ satisfy the equation 

^oir k,2n) = Ck - si, Ao(rj2„+i) = -Cfc - s^, /c G m, (1.26) 



they are real and the labeling is given by 



N 



<0 < ^fc,l < 4,1 < ^K2 < n,2 < -, k^ y , (1.27) 

4,0 < 4l ^ 4l < ^2 < 4,2 < 4i ^ 4,z-^ ^ " T' ^^'^^^ 

The resonances ?"^2n+s' s = 0, 1 (for A; G m, A; 7^ y) satisfy 

TTT? of 1 I IT 

4,n = (-^ ± ^fc,n)^ + go H , bk,2n = 0fc,O, &fc,2n+l = " " </>fc,l5 (1-29) 

0fe,s = I arccos — - — -^ — — G [0, |] as n — > oo. We describe the spectral properties of H. 

Theorem 1.5. (i) The periodic and anti-periodic eigenvalues A^ ^, fc G m = {1, .., m}, n ^ 
satisfy Eq. il.2^} and the relations M.2S\) - M.2f^) . 

(a) The resonances r^„, k Efn,n ^ satisfy Eq. M.2h]) and estimates \1.21^ - ^.2f^) . 
(Hi) For n ^ l,k Efn the following identities are fulfilled: 

aac{H) = Un-^iSn = U^aac{Hk), Sn = [^+_i, ^"J = lJT=0^k,n, n'^=0^k,n ^ 0, (1-30) 



^o,n = [A(l;„_i, XqJ, and A;.„ G (Tfe,„ = [r^„_ir;.„], (1.31) 

Gn = {E~,E+) = n^7A.,n, G2n = 70,2n. (1-32) 

(iv) If n ^ 1 is fixed, then G2n = iff In = ^- Moreover, \G2n\ = o(l) as n —>■ oo. 

(v) // y G Z, then 6*2™+! = /or all n ^ Q iff q & Le„en(0, 1). Moreover, each odd gap Gn 

has the form Gn = {,r~^, r+„) (n ^ 1 is odd) for some integer p = p{n) and r^„ satisfy 



vr^— + go ± kcnl + o(n ^) as n — *> 00, Qcn = / q{t) cosTmtdt. (1.33) 

4 Jo 



('wzj If ^ ^'L, then each gap G2n+i 7^ 0, ri ^ and |G2n+i| ^00 as n ^>- 00. 

Note that in the armchair case (see |BBKLj ) there exist non-real resonances for some 
specific potentials. That all resonances are real is a peculiarity of the high symmetry of a 
zigzag periodic graph. 

Finally, we formulate our results about the inverse problem. By Theorem 11.31 En = 

Hn^n ^ 1 are the eigenvalues of H and they satisfy (see [PT]) En = vr^n^ + Qo + >Cn-, where 
go = /o q(t)dt and (x„)f^ G £^ = i^. Here vr^ra^, n ^ 1 are the unperturbed eigenvalues 

and the real Hilbert spaces ip is given hy ip ={ f = ifn)T '■ Zln^i^-^^/n < oo>, p ^ 0. 

The monotonicity property Ei < E2 < ... gives that if g runs through L^(0, 1), then (x„)^ 
doesn't run through the whole space £^. In order to describe this situation, we introduce the 
open and convex set /C = {(x„)f G £^ : tt^ + Xi < (27r)^ + X2 < . . . } C £^. Let Hfn = Enfn for 
some nonzero eigenfunction /„ G 7i{En). Then /„^(^ 7^ for some u & Z and using Theorem 
11.21 we obtain 

/n,.(0) = /„,.(1) = 0, /;jO) = Gn,.., f'njlf = GI^^- > 0, /i„ G M, any n ^ 1, 

for some constant Gn,Lu 7^ 0. From Theorem II .3111 .41 we have a simple corollary. 

Corollary 1.6. (i) The operator H has only a finite number of non degenerate gaps Gn iff 
y G Z and q is an even finite gap potential for the operator —y" + qy on the real line. 

(a) The mapping $ : g ^ (go, (x„)^; {hn)f J is a real-analytic isomorphism between 

L2(0, 1) andRx IC x ij, where the set IC = {(x„)^ e i"^ : vr^ + xi < (27r)2 + x2< . . . } cf. 

(Hi) The mapping $e : 5' — ^ I ^o, i^n)T ] is a real-analytic isomorphism between Lg^g„(0, 1) 

and M X /C. 

For the convenience of the reader we briefly describe the structure of carbon nanotubes 
(see |Haj . |SDDj ) and explain how they are related to the graph F^. Graphene is a single 
2D layer of graphite forming a honeycomb lattice, as in Fig. ^ and |21 A carbon nanotube 
is a honeycomb lattice "rolled up" into a cylinder. In carbon nanotubes, the graphene sheet 




Figure 5: The unrolled honeycomb lattice of nanotube. The unit cell is spanned by the 
vectors ^i and ^2. The type of the nanotube is defined by the pair {Ni, N2) E N"^ , Ni ^ N2, 
and corresponding chiral vector f] = A^i^i + A^2^2- 

is "rolled up" in such a way that a so-called chiral vector ft = NiQi + A^2^2 becomes the 
circumference of the tube, where fli = 61 + 62, ^^2 = 60 + 61 see FigElandEl The chiral vector 
Q, which is usually denoted by the pair of integers (iVi, A^2), uniquely defines a particular 
tube. Tubes of type (A^, 0) are called zigzag tubes, see Fig and El They exhibit a zigzag 
pattern along the circumference, see Figure ^ (A^, A^)-tubes are called armchair tubes. 

Remark that in this zigzag standard physical model only each vertical edge has length 1 
and the other edges are shorter (or might be shorter). In this paper, we avoid this additional 
technical difficulty. Instead, we do our detailed spectral analysis for the simplest possible 
model of a physical nanotube T^ (to be described below and see Fig. EJ. This will simplify 
our analysis (in another paper |KLj ) of the Schrodinger operator J^b = {—i'V — £/y + q 
on T^ with a periodic potential q and a uniform magnetic field ^ = -6(0, 0, 1) G M^, 
B eM.. The corresponding vector potential is given by £/{x) = |[^, x] = ^{—X2,Xi, 0), x = 
{xi,X2,X3) G M.^. Maybe the second simplest model is the armchair nanotube [Ha], but then 
the monodromy matrix is 4 x 4 and there exist complex resonances for some specific potential 
q (see [BBKL]). 

Finally we consider the Schrodinger operator J^ = —9^ + q with a periodic potential q 
on the zigzag nanotube T^ C M^. Our model nanotube T^ is a union of edges T^ of length 
1, i.e., T'^ = UujfzzTu}, see Fig. HI Each edge T^ = {x = r^ + te^, t G [0, 1]} is oriented by 

We have the coordinate x = r^ + te^ 

n,j,k) e Zhy 



tu 



the vector e^ G M and has starting point r^^ G 1 

and the local coordinate t G [0, 1] (length preserving). We define r^, e^^, 

3n 

^n,0,k = ^n+2fc + ~^^05 ^n.l.fc = T^n,0,k + Gq, T^n,2,k = 1^n+l,0,k, 

_ _ , *^0 _ ^0 
^n,0,k — GQ) Gn,l,k — >^n+2k+l ~ ^n+2k + "i^j Gn,2,k — »^n+2k+2 — ^n+2k+l ^, 
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where x,- = Rn{cj,Sj,0) G M^, -Rat = . _;^ , The points ro,o,fc are vertices of the regular 

2N 

N-gon <^o- The vertical edge To,o,a: is the segment and lies on the cylinder C = {(xi, X2, x^) G 
M^ : xf + x^ = Rn}- The starting points ri_o,A: = i"o,2,fc = ><^i+2fc + |eo, k G Z^ are the vertices 
of the regular N-gon ^i. 3^\ arises from ^o by the following motion: rotate around the 
axis of the cylinder C by the angle j^ and translate by |eo. The non-vertical edges T^x^ 
and To,2,fc have positive and negative projections on the vector Bq. Repeating this procedure 
we obtain all edges of T^ . Note that each non- vertical edge TQjk,j = 1,2 (without the 
endpoints) hes inside the cyhnder C. 

Our operator J^ on T^ acts in the space L^(T'^) = (BuiezL^iT^)- Then acting on the edge 
T(^, J^ is the ordinary differential operator given by {J^f)uj = — /" + qfuj, fuj, f" G L'^{T^), 
to e Z where q G L^{0, 1) and / G Tl{J^) = W^{T^) satisfies the Kirchhoff Conditions 
fll.l|) . Thus the operator J^ on T^ is unitarily equivalent to the operator H on F^. 

We remark that such models were introduce by Pauling paj in 1936 to simulate aromatic 
molecules. They were described in more detail by Ruedenberg and Scherr [RS] in 1953. For 
further physical models see [Ha], [SDD]. For papers on spectral analysis of such operators, 
see the references in JKuj. 

We present the plan of the paper. In Sect. 2 we construct the fundamental solutions 
and describe the basic properties of the monodromy operator. Moreover, we prove Theorem 
1.2 about the eigenf unctions. In Sect. 3 we prove the basic properties of the Lyapunov 
function Aq which is entire on the plane. This function is important to study the functions 
Afc, k G Zat which are analytic on the two-sheeted Riemann surface TZk- In Sect. 4 we prove 
the basic properties of the Lyapunov function A^, k G Z^r and Corollarv 11.61 

2 Fundamental solutions and eigenfunctions 

Proof of Theorem 11.11 (i) Define the operator S in C^ by Su = {un,ui, . . . ,un-i)^ , 

u = {un)i G C^. The unitary operator S has the form S = "^^ s'^Vk, where Sv^ = s^Vk 
and Vk = -^(lyS'^jS'"^^, ...,s~''^~^^) is an eigenvector (recall s = e*"^); VkU = Vk(u,Vk) is 
a projector. The function / in the Kirchhoff boundary conditions (jl.ll) is a vector function 
/ = ifuj),^ = {n,j,k) G Z. We define a new vector- valued function fnj = {fn,j,k)k^ 
where each /„,,•, {n,j) E Zi = Z x {0, 1,2} is an C^— vector, which satisfies the equation 
~fn,j + lfn,j = '^fnj, and the conditions (which follow from the Kirchhoff conditions (jl.l|) ) 

/n,0(l) = /n,l(0) = 5/„,2(l), /„+l,o(0) = /„,l(l) = fn^, (2.1) 

-/;o(i) + /;i(o) - sji,{i) = 0, fUi,oio) - /;i(i) + 42(0) = 0, (2.2) 

for all 72 G Z. The operators S and H commute, then H = ®i {HVk}- Using ()2.H) . 
fl2.2|) we deduce that HVk is unitarily equivalent to the operator Hk- The operator Hk on 
F^ = UaezJ^\ acts in the Hilbert space L^(F^). Then acting on the edge F^, the operator 
Hk given by: {Hkf)a = -fa + (l(t)f^, where /„, f'^ G ^^(0, 1), a = {n,j) G Zi, on the vector 
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functions / = {fa)aGZi, which satisfy the boundary conditions p.4|) . (jl.Sp . i.e., 

/n,0(l) = /m(0) = sVn,2(l), /n+l,0(0) = /n,l(l) = /n,2(0), S = e^^, (2.3) 

-/;o(i) + /;i(o) - s'ji,{i) = 0, /;^,,o(o) - /;i(i) + /^^(o) = o. (2.4) 

(ii) Fix k G Z^. We determine the fundamental solutions 9fe(x, A), $fc(x, A) of the equa- 
tion -/" + qf = \f, A e C, on T\ where 6^, $fc satisfy (|T31)-(j23I) and 

efc,„(0, A) = <l>',,„(0, A) = 1, e'.jO, A) = $fc,,(0, A) = O, a = (0, 0). (2.5) 

Below we assume / = 6fc or / = <l>fc and let {It = '(9(1, A), ^9^ = ip(t, X), .... Any solution of 
the equation -/" + qf = \f satisfies /(t) = 79t/(0) + f^(/(l) - ^i/(0)), t G [0, 1] and 

Substituting (j2.6j) for / = /oi, / = /02 into the first Eq. in (j2.4|) at ra = 0, we obtain 

-^1/0,0(1) + (/o,i(l) - ^1/0,1(0)) - ^'=(^'1/0,2(1) - /o,2(0)) = 0. 
Then the condition (j2.3p gives 

<^i/o,o(l) + (/i,o(0) - ^1/0,0(1)) - s' {^[s-'f o,o{l) - /i,o(0)) = 0, (2.7) 

which implies 

-¥.1/^0(1) + m/i,o(0) - 2A/o,o(l) = 0, M = 1 + s^ = 2cksl (2.8) 

Thus we obtain the first basic identity 

/i,o(0) = ^/;o(l) + — /o,o(l), u = 2cksK Cfc = cos^. (2.9) 

Substituting ()2.6p for / = /oi, / = /02 into the second Eq. in ()2.4|) at ri = 0, we get 

^ifioW - iv'JoAl) - /o,i(0)) + (/o,2(l) - ^1/0,2(0)) = 0. 
Using ()2.3p . we obtain 

VifloW - (¥^'1/1,0(0) - /o,o(l)) + (s-Vo,o(l) - ^1/1,0(0)) = 0, (2.10) 

which yields 

V?i/{.o(0) = 2A/i,o(0) - ^2/0,0(1), U2 = l + s-'' = 2ckS-l (2.11) 
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Then substituting ()2.9|) into ()2.1ip . we obtain the second basic identity 

/(,o(0) = ^^^' ~ "' V o,o(l) + — /o,o(l), n = 2c,si. (2.12) 

ULpi U 

Substituting consequently f = Qk and / = $fc into ()2.9p and into ()2.12p . and after this 
0fc,.(O, A), e;, JO, A), .. into MkiX) = ("gf' jjj; J| |f' Jq^ J)) , « = (1, 0), we obtain (d. 



Using (fOl) we obtain det 7^ = ^-r^ = s '', det Mk = det T^AI = s '', since det A^ = 1. 
Moreover, using TZAiTZ^^ = ( „/ / I , ^ = I ^ I , we obtain 



Tr A<. = liTkTZMTZ-' = ^-^ Tr 



_i_2._2 / 2A iW^i 1 

c, Ua2 - Aci 2A i Ui^; y,; 



2^ 



= 2A^i+y.i#i + 4 A2-4 +2A^; = 2A2 + ^-cn = ^^ 

Cfc V / Cfc V 4 y Cfe 

which gives (jl.lUp . ■ 

Recall that if°,S)(/7°) = W\T^) is self-adjoint [Cal]. Let \\ff = /^^ |/(x)pda; and let 
ll/^llo = /o^ \fUt)?dt for /^ G L2(r^),cj G Z. Assume that 

\\qff^^\\H'ff + C\\fr, all fe^iH'), (2.13) 

for some constant C > 0. Then by the Kato-Rellich Theorem (162 p. [ReSl]), H = H^ + q 
is self-adjoint on S}(if'') and essentially self-adjoint on any core of H^. We prove ()2.13|) . 

Lemma 2.1. Let q G -^^(0, 1). Then the estimate \2.1'J\) holds true. 

Proof. For / G S)(if°), / = {foj)ujez,^ = {n,j, k) G 2, we have the following estimates 

||g/J|o^max|^(x)|||g||o, max|^(x)|^ / (|/.| + |/:i)rfx ^ l|/:i|o + l|/J|o, 



which yield 

Thus for Cq = \\q\\l we obtain 

\\qfr ^ 2c,j2(\\ai + ii/jio) = 2CM'r + ii/ir) ^ ^ii^vir+(^ + 2c,) 

for any e > 0, since 

2C,||/'f = 2Cq{H^f, f) ^ 2Cq\\H^f\\ ll/IK eWH'^ir + ^1' "'" 
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which yields ()2.1|) . ■ 

Note that the operator Hk in L'^iV^) with S)(/7fc) = {/, ®a&zj'a e ^^(r^) : / satisfies 
the Kirchhoff conditions ()1.4j) . (jl.5j) } is self-adjoint. The proof is similar to the case H. 
Proof of Theorem 11.21 (i) Let HkF = XF for some eigenfunction F and some A G ctd- For 
eachF,^ = F\r^,uj = {n,j) G Zi = Zx J, we have F<^(t) = F^{0)+J^ F^{s)ds,t e [0,1]. Then 

|-F1j(0)| ^ |Fa;(x)| + Jg \F^(t)\dt and integrating over x G [0, 1], using the Holder inequality, 
we obtain 

\F^m< [ i\FUt)\ + \Kit)\)dt^\\F4o+\\Fjo = oil), as n ^ ±oo, (2.14) 

since F G 7ifc(A). Furthermore, each restriction F^ has the form F^{t) = a^j^t+buj^Pt, t G [0, 1] 
for some constants a^;, b^, which implies F^{1) = F^(0)'i9i, where ipt = (p{t, X),^t = ^{t, A), .. 
The Kirchhoff conditions give F,+i,o(0) = i^„,i(l) = ^iF,,i(0) = t9iF,,o(l) = ^?i^n,o(0), 
which yields 

Fo,o = t9f F_„,o(0), all n e Z. (2.15) 

Let |i?i| ^ 1, the proof for the case l-i^il > 1 is similar. Then ()2.15|1 . ()2.14|1 imply Fq^ = 
i??"F_„,o(0) = o(l) as n ^ +oo. Thus flTT^ gives F„,o(l) = for all n e Z. Finally, F 
vanishes on all vertexes of F^, since the set of all ends of vertexes F„_0; n ^ Z coincides with 
the vertex set of F^. 

(ii) Using (ITTTIl - drT^ . we deduce that ip^^'^ satisfies the Kirchhoff conditions dHj), (fT3|) . 
Thus -0° is a eigenfunction of H^. 

The operator H^ is periodic, then each tp^"'\n G Z is an eigenfunction. We will show 
that the sequence ■ip^"'\n G Z forms a basis for 7ik{\). The functions ^^"-' are linearly 
independent, since supp ■?/'*•"''' H supp ■i/'*^'"^ = for all n 7^ m. 

Consider the first case rj = 1 — s^ip\ 7^ 0. For any / G 'Hk{X) we will show the 

identity dni, i.e., / = / , where / = E„ez /n^^"^ L = ^- The definition of /and 

A G ctd give /|r„o = /lr„o = /"V^ fo^' all n G Z. This yields Y. |/„p < 00 and / G -^^(F^), 
since / G F^ij:^) 

Note that / satisfies the Kirchhoff conditions (jl.4j) . (jl.5j) and — /^' + qf^ = Xfa, « G Zi. 
Consider the function u = f — f = {ua)a£Zi, where Unj = Cnj'^it, X),j = 1, 2 and Unfl = 
for some constant Cn,j,n G Z. The Kirchhoff boundary conditions ()1.4|l - ()1.5j) yields u = 0. 

Consider the second case s''(f'i = 1. For any / G ?ifc(A) we will show the identity 
(HHSl), i.e., / = /, where / = Enez/"^^"'*' /" = /n,i(0)- ^^°"^ ^^^ definition of / and 
A G (Td we deduce that /|r„,i = /|r„,i = fnf for all n G Z. This yields ^ |/„p < 00 and 
feL\r^), since / G ^^(F^). 

Consider the function u = f — f. The function u = at all vertices of F^ and then 
Ua = Ca^Pt, C( = {n,j),n G Z, j = 0, 2. Assume that Cq^ = C. Then the Kirchhoff boundary 
conditions p.5|) yields Cno = —Cn 2 = —C and C„ = Cn+i = C, which give C = since 
mGL2(F«). ■ ' ' 
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3 The Lyapunov function Aq 

Recall that 'd{x, A), ip{x, A), a; G M are the fundamental solutions of Eq. —y" + qy = Xy, A G 
C, on the real line such that ?9(0, A) = ip'{0, A) = 1, i9'{0, A) = (^(0, A) = 1. For each x eR 
the functions "d, •&', ip, if' are entire in A G C and satisfy: 

i^{x,\) =cos^/\x + ^^ I sin \^x + sm^/X{x- 2t)]q{t)dt + 0l — -- — J, (3.1) 



/(x,A) = -\/Asinv^x + - / UosV\x + cos^/\{x-2t)jq{t)dt + 0^-—J^], (3.2) 

Sinv/Ax IT/' r- r- \ /g|ImVA|a;\ 

(^(a;,A) = ^~ + ^/ (-cosVAx + cosVA(a;-2t)jg(t)dt + Of — ^y^j, (3.3) 

1 /"^ / \ / ^\lm'/\\x\ 

ip'{x, A) = cos yXx + ^^ / I sin VAx - sin v^(x - 2t) ) g(t)rft + O — -- — (3.4) 

2vA Jo ^ ^ V 1^1 / 

as |A| ^ cxD, uniformly on bounded sets of (x; q) G [0, 1] x L^(0, 1) (see [PT]), and 

r- gosinv^ 0(eli^^l) /"^ . ^ , 

A(A) = cosv^+ ^° ^^^ ^' ^"^y ^^^^ ^ ^ ^ 

Substituting dSHJ-dSini) into Aq we obtain (frT9|l . i.e., 

0,-, 9gosin2yA /e^li'^^v^lA 9cos2v^-l 



A„(A) = AS(A) + ^^^ + O^^J, Ag(A) 

Let A_ = |(/(1, •) - ^(1, •))• Substituting the identity </?(!, ^^'(l, ■) = A^ - A^ - 1 into 
Ao = 2A2(A) + ^'^^-y^'^) - 1, we obtain 

Ao^^^^^^, where ^_ ^Ui_^M. (3.6) 

Asymptotics ()3.1|l - ()3.4j) yield as |A| ^00 

A_(A) = ^ + -^^- ^, F{X) = sinVXil- 2t)q{t)dt. (3.7) 

2vA |A| Jo 

Proof of Theorem 11.31 (i) Using ()1.10j) we obtain the characteristic equation 

det{Mk-Tl2) = T^ -2akT + s-'' = 0, au = — ^ = ^-^, ^^ = Aq + 4- 

2 1 + s*^ 
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The eigenvalues Tk^± oi Mk are given by Tk,± = a^ ± a/o^ — s^^. Using ak ^^ ' = ^k and 
^l ~ ^'^) — Pk we have 



0-f!nif^2 ,-fc 



1. . -1. 1. . . . n + s'^) (1-s'^' 



Afc,+ = -(rfc,+ + r, ^) = -(rfe,+ + sVfc,_) = a^ ^ ^ + ' \ Ul -s-^ = ^k + ^k 



and a similar argument yields ^k- = ^k — y/Pk- 

(ii) The standard arguments (see [Cal]) yields a{Hk) = aac{Hk)Uapp{Hk), i.e., asc{Hk) = 
and aac{Hk) = {A G M : Afc(A) G [-1, 1]}. 

(iii) Let A'^(Ao) = and Afc(Ao) G (—1,1) for some Aq G cTk,k = 0,..,N. Then we 

have the Tailor series Afc(A) = Afc(Ao) + C^^^^^ + 0{Cp+^) as C = A - Aq ^ 0, where 

Ajf (Aq) 7^ for some p ^ 1. By the Implicit Function Theorem, there exists some curve 
F C {A : |A-Ao| < e}nC+,F 7^0, for some £ > such that Afc(A) G (-1,1) for any A G Y. 
Thus we have a contradiction with ()1.17|) . ■ 
In order to prove Theorem 11.41 we need 

Lemma 3.1. Let q G L'^{0, 1). Then the following statements hold: 

(i) Ao(A) ^ 1 for all X = fin and X = z/„_i,n ^1. If in addition A'^{X) = 1, then Ao(A) = 1. 

(zz) IfA{r]) = for some t] eR, then Ao{r]) = -^±fShnl ^ _|. 

(zzz) Let q G Ll,^^{0, 1) = {g G Ll^^^{0, 1) : q{x)=q{l-x),xe[0, 1]}. Then 

QA^ — ^ 
Ao = ^^, A = y,'(l,-)=^(1,-), (3.8) 



7o,2n = 7n, where each 7„ = (/i„, z/„), or 7„ = (z/„,/i„), n ^ 1. (3.9) 

r/ie zeros ?7„ = Ao,2n-i for all n ^ 1. 

(iv) There exists an integer rio > 1 such that Aq — 1 has exactly 2no + 1 roots, counted with 
multiplicities, in the domain {A : \\/X\ < 7r(no + |)} and for each n > uq, exactly two roots, 
counted with multiplicities, in the domain {A : |vA — nn] < j}. There are no other roots. 
(v) There exists an integer uq > 1 such that Aq + | has exactly 2nQ roots, counted with 
multiplicities, in the domain {z : \\^\ < imo} and for each n > no, exactly one simple root 
in the domain {A : |vA — 7r(n — |)| < ^}. There are no other roots. 

(vi) There exists an integer uq > 1 such that Aq has exactly no roots, counted with multi- 
plicities, in the domain {A : |vA| < tcuq} and for each n > Uq, exactly one simple root in 
the domain {A : | vA — 7r(n — |)| < ^}. There are no other roots. 

(vii) Let c G (— |, 1) and let u^ = nn ± Uq, Uq = ^ ^ G [0, |], n ^ 1 and Uq < u^ < 

Ul < U2 < ■■. Then there exists an integer uq > 1 such that Aq — c has exactly 2no + 1 roots, 
counted with multiplicities, in the domain {X : |vA| < R},R = {u^^ +u~^^i)/2 and for each 
n > no, exactly one simple root in the domain {A : | vA — u^\ < r for some small r G (0, 1). 
There are no other roots. 



Proof (i) Let A G {/i„, u^} for some n ^ 1. Then we have Ao(A) = 2A2(A) + ^'^^■^^'^^^'^^ - 1 
2A2(A) -1^1. Moreover, if A G {/i„, z/„} and A2(A) = 1, then we have Ao(A) = 1. 
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(ii) We have = 2A(r7) = (/?'(!, r/) +-(9(1, r^), which yields ^^'(l,?]) = —{}{l,ri). Substituting 
the last identity into Ao(A) = 2A^(A) + ^ — 1 we get the needed estimate. 

(iii) It is well known that if g e i^Ln(0' 1)^ then A = (^'(1, ■) = ^(1, ■), see p.8, [MW]. It 
is well known that if A = ^'(1, ■) = ^(1, ■), then q G Ll^eni^, 1), see [PT]. 

(iv) Define the contours C„(r) = {A G C : |vA — 7rn| = vrr}. Fix another integer rii > uq. 
Consider the contours Cq^uq + (1/2)), Cq^ui + (1/2))), C„(l/4), n > tiq. Using the estimate 
g|im\/A| ^ ^1 gjj^ a/A|, I a — 7m\ ^ ^ and (jl.lQj) we obtain on all contours (for large no) 

|(Ao(A) - 1) - (A°(A) - 1)1 = 0(6^11- -^1) = I sinv^lVl) = o(l)|A°(A) - 1| 

Hence, by Rouche's Theorem, Ao(A) — 1 has as many roots, counted with multiplicities, as 
Aq(A) — 1 in each of the bounded domains and the remaining unbounded domain. Since 
Aq(A) — 1 has exactly the simple root = and one root of the multiplicity 2 at each vr^n^ ^ 1, 
and since rii > uq can be chosen arbitrarily large, the point (iii) follows. 
The proof of (v) and (vi) is similar. ■ 

Lemma 3.2. Let c G [— |, 1]. Then the equation Ao(A) = c, A G C has only real zeros, which 

satisfy 

(i) If c & ( — |, 1) and go = /o q(t)dt, then these zeros z^ are given by 

± ± , ?o , o(l) 



zt < z^ < z'l < Z2 < .. and y z^ = vr + —-— -\ — as n — ;> 00, (3.10) 

2ut n^ 



arccos "^ 



where u^ = im ±: u^, n ^ 1 and u^ = ^ ^ G [0, |], Uq < u^ < u'l < U2 < ... 

Moreover, the zeros z^ have another forms given by: x~ = z^_^_i,x^ = z~, and 

Xi < x^ < X2 < X2 < x^ < .. \f^ = "^n + TTT H T' ^^ n — i> 00, (3.11) 



2v^ n 

where v^ = n:{n - 1) ± (| - u'^). 

(a) If c = 1 or c = —^ . Then these zeros z^ and x^ are given by 

Zn < z7 ^ zi < z^ ^ zf < .. and \/z^ = wn -\ 1 —, if c = 1, 

_ , _ , /-^r Ix Q'o 0(1) 5 

x^ ^xj < X2 ^xj < .. and ^^x^ = n{n - -) + \ —, if c = --, 

2 27m n^ 4 

as n ^ 00. 

Proof, (i) Let z = vA. Recall the folowing asymptotics from [KKl]: 

A(A) = cosx(A), x(A) = v^-^ + ^ as |A| ^ 00. (3.12) 

vq^i^m qa2 a2 C 

Using ()3.6p we rewrite the Eq. Ao(A) = c, A G M in the form c = -^ — , which yields 

cos2x(A) = A, + ^^y^, A, = i±^G(-l,l). (3.13) 
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In the case g = we have cos2z = A^ and the corresponding solutions \/z^ = u^. Then by 
Lemma ITTl \^/z^ — u^\ < r for all n > uq, for some uq ^ 1 and some small r G (0, 1). Using 
fl3.13|) . ()3.7|) we obtain cos2x(A) = Ac + o{X~^) as A -^ oo, and then (i^ = \f^ — m^ — > 0. 
Thus e^ = C^ — -^ + ^V- -^ and the Taylor series gives 

cos2x(A) = cos2(M^ + e^) = cos2m^ - sin(2M^)£^(l + 0(e^)) = Ac + o(n~^). 

Thus we obtain £± = (±--2^ + ^ = ^ which implies, C± = -*^ + ^. 

Proof of (J3.1H) and (ii) is similar. ■ 

Proof of Theorem 11.41 (i) The function Aq = 2A^ H ^ '''^^ "'' — 1 is entire and real 

on real line. Moreover, it has asymptotics (J1.19p . By Theorem 11.31 the function Aq has only 
simple real zeros riQ^n,n ^ 1, which satisfy ?7o,i < riQ2 < ''7o,3 < •••• Then by the Laguerre 
Theorem (see Sect. 8.52 [Ti]), the function Aq has only real simple zeros Ao^n;'"- ^ 1) which 
are separated by the zeros of Aq: ?7o,i < Ao,i < ?7o,2 < Ao,2 < Vo,3 < ■■•; since Ao(A) ^ oo as 
A -^ — oo. Moreover, we obtain (— l)"Ao(Ao,n) ^ l,n ^ 1. 

Define an interval Q = (— oo, tt^tIq). Lemma f3 . II gives for large integer uq > 1 : 

A) The function Aq — 1 has 2no + 1 zeros on Q. 

B) Ao(r7o,n) ^ -|, ^ = 1, .., no, where 77^ G ^ is a zero of A, i]i < 'r]2 < .... 

C) Ao(/i„) ^ l,n = 1, ..,no, where fin ^ Q are zeros of ip{l, A), the Dirichlet eigenvalues. 
These facts and fii < rji < fi2 < ri2 < ... yield /i„ G 7o,2n and rj^ G 7o,2n-i for all n ^ 1, 

where 7o,n = [A(^„,A^„]. Moreover, we obtain (jl.2(jp . 

(ii) Using (i) we deduce that the periodic and anti-periodic eigenvalues A^„,n ^ satisfy 
fll.l8|) . Consider A^2n+i5 which satisfies Ao(A) = —1. Using Lemma |S!21 (i) for c = — 1 we 
have the first asymptotics in p.2H) . In order to show p.2H) we consider A^2n with even 
n ^ 1. The proof for odd n is similar. Using ()3.H) - ()3.4|1 and ()3.7|1 we obtain 

Oil) - 1 + ^^ 0+^^. 0+^^ 

Ao(A) = ^, Ao(A) = -^±^, A_(A) = ^^^^^, A_(A) = -^^^^^, (3.14) 

0(1) - 1 + ^^ ... 0(1) • dA 

A(A) = ^, A(A) = -^-^, A(A) = ^, where A = — (3.15) 

as vA = Tm + 0{l/n), where Qcn = /q (1 ~ 2t)g(t) cos2TTntdt. Using Ao(Ao,2n) = 0, we have 
Ao(A„) = Ao(Ao,2n)s„(l + 0(s„)), s„ = A„ - Ao,2n as n ^ 00. 

Q A 2 A 2 cr 

The identity A(A„) = and Aq = f- and (ITTHl yield 

. A_(A„)A_(A„) _ (g.n + 0(^))(gcn + 0(^)) .„ ^ „. 

Ao(A„) j^^, . (3.16) 

Thus asymptotics ()3.14|) - ()3.16p give 

^ _. . _ fen + 0(^))(9cn + Q(^)) 

Sim 
18 



which imphes that roughly speaking the point Ao,2n is in the center of the gap 7„, since 
A„ = (7m)2 + go + 0{l/n) see [Ko]. 

We will determine asymptotics of 7o,2n = (-^o 2n' ^t2n)- -Due to (j3.15|) we get 

A(Aj2J = A(A„)+Aj„, ^t = A(A„)^(l + ^^), Co':2n = AJ2. - An, (3.18) 



1 = A(A^)=A(A„) + A^, A^=A(A„)^(1 + ^^), e = A^-A„, (3.19) 



which yields A(AJ„) - 1 = -A^ + AJ„. Substituting ^^^T^ into the identity A(AJ„) 
1 + ^^^^ - 1 we obtain 



^^^°-) - ^ - V ^ ^ " 18(2vrn)^ ' ^^'^^^ 

Substituting (jT^ and (Hnil into the identity A(Aj„) - 1 = -A^ + Aj„ we get 

e 0(|g„|) _ ±2 ^ 2 0(1) 

Using A^ = (7m) ^ + go ± l^nl + 0{n~^) from [MO] and A„ = (7m) ^ + go + 0(n"^) from [Ko] 
we obtain C^ = A^ — A„ = ±|g„| + 0(r2~^) which yields 



Cot = ± V^" ' - f + ^ = ^V '^'^'^ - f + ^^""^- ^'-''^ 

(iii) Let /i„, z/„ lay on the different edge of the gap 7„. Then Lemma lXTl fi) gives 7„ = 7o,2n- 
Conversely, let 7„ = 7o,2n = (o, &) and let A G {a, 6}. Then the identity Aq = 2A^ + 

tl^^iMhl _ 1 yields 1 = Ao'(A) = 1 + ^'^y^\ thus ?9'(A)<^(A) = 0, which gives (iii) 

(iv) Let 7„ = (a, 6) and let a < /i„ < i^^ < 6. The proof of other cases is similar. Let 

Ao = /o + /, where /o = 2A2 - 1, / = ^'^^'y^'^) . We obtain 

/ol7n > 1> /U > 0, /l7„v < where a = (/i„, z/„), 

which yields AqIo- > 1 and Ao(A) < 1 for any A G {a,b}. Then there exist two points 

A^ G (a, 6) such that Ao(A^) = 1. Thus we have 7o,2n C 7„ for all n ^ 1. 

We will show 7o,2n = iff 7n = 0- If In = 0, then 7o,2n C. 7„ for all n ^ 1 yields 7o,2n = 0- 
If 7o,2n = 0, then for A = {/i„, z/„},n ^ 1 we obtain Ao(A) = 2A2(A) + ^'^^^^'>f^^^^ - 1 = 

2A2(A) -1 = 0, which yields A'^{X) = 1. Thus we deduce that 70,2™ = iff 7„ = 0. 

(v) Let g G L^^,„(0, 1). Then we get ip'{l, ■) = ^(1, ■) = A (see p. 8, [23]), which together 

with dSini) yields Aq = f A^ - |. Then the zeros of A(A) and Ao(A)(at Ao(A) < 0) coincide, 

since all zeros of Aq are simple. 
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Conversely, let Ao(Ao,n) = — | for all odd n ^ 1. Then Lemma EH] ii) gives -(9(1, Ao,n) = 
^(Ao,n) = '^'{^,^o,n) = for all odd n ^ 1. Then rjn = '7o,2n-i5''^ ^ 1, which implies 
(p'{l, ■) = 'i?(l, ■) = A. Thus the Wronskian identity gives A^ — 'i9'(l, ■)v^(l, ■) = 1 ^^^ all 
zeros of A^(A) = 1 coincide with the Dirichlet and Neumann eigenvalues. Then the results 
of [GT] or [KK] imply q G Lg^g„(0, 1). It is simple fact in the inverse spectral theory and it 
can be proved using other methods, see [PT]. ■ 

4 The Lyapunov function A^^, /c = 1, .., m 

Proof of Theorem 11.51 (i) We determine the equation for periodic eigenvalues. Using 

()1.10|) and 44 = (1 + ■s'')(l + s~^) ior k em = {1, 2, .., m}, we have 

k , , ^-k^Ao + 4 4(c2t(Ao + 4)) 



det(A^fc ^ /2) = 1 ^ TrMfc + s-*= = 1 + 5-^= ^ 4 



which yields the equation 

2iTk 
Ao(A) = cl — si = cos — - G (—1, 1) \ {0} for periodic eigenvalues A, (4.1) 

Ao(A) = —cl — si = —1 for anti — periodic eigenvalues A. (4.2) 

Then by Lemma 13.21 all zeros of det(A^fc =F 12) = are real and simple. Using ^^ = 

2 

Ao + si, pk = ^{cl- Q, we get (recall ii = j^u) 

2v^ V cl ^^ 
For periodic eigenvalues A = A^„, where n ^ is even, we have Afc(A) = 1 and 

ik = Ao + sl = cl-sl + sl = cl, v^ = I - ^f^ = I - d = si, at X = A^„, (4.4) 
and then ()4.3|) yields at A = A^^: 

^'-^ ^^^^-1-1=0, A..=-^rA-i^u-^ri+4^<o, 



Afc = a + v^ = Ao + 4 + v^, Afc = a + ^^ =( 1 - :|^ ) Ao. (4.3) 



Ao 4v^ ' ' c| VVPfc 2pl^V 4\ 4 

which gives ()1.23|) for periodic eigenvalues X^2n- 

For anti-periodic eigenvalues A = A^^, where n ^ 1 is odd, we nave Afc(A) = —1 and 

^„ = Ao + 4 = -l + 4 = -cl, ^=-l-^, = -sl, atA = Aj„, (4.5) 

and then ()4.3|) yields 



Ao '44 ' 4 \Vp^ 2pl^y 4 V *■ 



k 
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which yields ()1.24p for anti-periodic eigenvalues A^^. 

Let A = A^2n5 ^ G m, n ^ 0, be periodic eigenvalues. Recall that they satisfy Ao(A) = 
cos ^ G (-1, 1). Using Lemma K^ for c = cos ^ we have (fT^ . 

(ii) We determine the equation for resonances for k Gfn. We have: Pa;(A) = <;^ ^fc(A) = 
±Cfc. Consider the first case +: n is even, i.e., ^ G Z. We obtain 

6(A) = Cfc ^ Ao(A) >a+ = Cfe -4 G (-1,1), k Em. (4.6) 

Then the resonances r^„ are zeros of Eq. Ao(A) = Ck — si E (—1, 1), A G C From Lemma 
111 21 we deduce that all these resonances are real. Moreover, we have 1 > Oj*" > a^ > ... > a^, 
which yields for even all even n ^ 1: 



^1,0 ^ ''^2,0 ^ ••• ^ ''"m,0' ^m,n--- < ''"2,n < ''"l,n < ^q^^ < Aq„ < T-,^,^ < r2„ < ... < 



r' 



7m,n C 7m-l,n C ... C 7l,„ C 7o,n, Ct^C? Gn = n^7fc,„ = 7o,n. (4.7) 

Consider the second case — : n is odd. We have 

5 A^ 

6(A) = -Cfe ^ Ao(A)=a^ = -Cfc-4G(--,-l), j ^ Z, (4.8) 

5 A^ 

efc(A) = -Cfc ^ Ao(A) = a^ = -Cfc -4 G [--,-1), y G Z. (4.9) 

Then from Theorem 11.31 we deduce that the zeros of pk are zeros of Eq. Ao(A) = — c^ — si E 
[— |, —1), A G C. Moreover, by Lemma f3. 21 they are real and simple. 

If y ^ Z, then y = p + e for some integer p ^ 1 and e G (0, 1). Then we obtain 

70,n 3 7l,n 3 72,n 3 ••• 3 7p,n, and 7p+l,n C 7p+2,n C .. C 7m,n, (4-10) 

7p,n C 7p+i,n or 7p,„ D 7p+i,„ anrf G„ = r\^jk,n = 1p,n H 7p+i,n, (4-11) 

for all odd n ^ 1. If p = y G Z, then we get for all odd n ^ 1: 

70,n D 7l,n ^ 72,n D ••• D 7p,n, Onrf 7p,„ C 7p+l,n C .. C 7m,n, (4.12) 

Gn = n-=o7M = 7p,n, Ao(rJJ = -^. (4.13) 

Resonances ri^2n,kEfn,n^0 satisfies Ao(r^2n) = c^ — s^ G (—1, 1). Using Lemma IT^ for 
c = Cfc — s^ we obtain asymptotics ()1.29|) for r^2n- 

Resonances ^^2n+i'^ G m. A; 7^ y,n ^ satisfies Ao(r^2n+i) = ~"Cfc - s^ G (-|,-1). 
Using Lemma IH^ for c = c^ — s^ we obtain asymptotics ()1.29p for r^2n+i- 

(iii) Finally, we obtain the following identities for (n, k) eN xm: 

Sn = U™^o'^fe,n, O-0,n = (Ao,„-l, Aq^^), (Tfc,„ = ('^fc.n-l' '^fc.n)' (4-14) 
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Gn = n^LoTfc.n, 70,n = (Ao,„, Xt,n) ^ lk,n = K,„,'^J„), (4.15) 

which together with (jHH), (jUinil-dlllSl) yield gives (fT3n|l - ffn!^ . 

(iv) Using Theorem II .41 we deduce that G2n = 7o,2n = iff 7n = 0. Moreover, asymptotics 
|7rt| — > as n -^ cxD give \G2n\ ^ as n —^ oo. 

(v) In the case p = y G Z and odd n ^ 1 using the identity ()4.13|) . we have Gn = jp^n = 
(r~„,r+„), where r^„ are zeros of Eq. Ao(A) = — |. Theorem 11.41 (v) gives r^„ = A„ iff 
gGLL„(0,l). 

We determine the asymptotics p.33|) . Let (z^)^ = r^n- Using ^ — = ^ we have 

9A^(r^„) = A?.(r^„). Thus in the case g = we have cosz = and the corresponding zeros 
are given by Zn = vr(n — |), n ^ 1. 

In the case g 7^ we have A(A) = ±|A„(A) = %^ as A ^ 00. Then by Lemma l3.H 
\z^ — Zn\ < J for all n > uq, for some uq ^ 1. Thus we obtain z^ = Zn + Cn^Cn ~^ 0- 
Moreover, using ()3.12|) . we get ^ = A(z„ + (^) = cos(z„ + Cn ~ ^ + ^)' which implies 
Cn = ^ + ^- Let C^ = 3^,vt^0 and for i?± = i| + ^ using ^^ we obtain 

A(^„ + O = cos(z„ + R^) = i-ir-' cos(-^ + R^J = (-1)"-^ sin(^ + ^^^^ ^ 



2 «/ ^ / V22„ n 



2 



Using this and ()3.7p we get 



gc^^o(i) 

22„ n^ 



A-(rJJ = (-ir^ + ^, ^-00. (4.16) 



Thus the Eq. A(r± ) = ±iA_(r± ) yields (IT3!^ . 



(vi) In the case y = p+e, where p ^ 1 is integer and e E (0, 1) and odd n ^ 1 the identity 



TV 



(j4.11|) we have (?„ = 7p_„ fl 7p+i,n- In (!1.29|) we will show that |7fc,n| ^ooasr7,^oo,A;7^ ^ , 
which yields IG^I ^ 00 as n ^ 00. ■ 

We will prove the last Corollarv 11.61 

Proof of Corollarv 11.61 (i) Recall that by Theorem 11.51 each G2„ = 7o,2n C 7„,n ^ 1. 
Thus by Theorem 11.51 fiv). G2n = 7o,2n = for all n ^ uq and for some tiq ^ 1 iff g is a finite 
gap potential for the operator —y" + qy on the real line. 

Consider the odd gaps G2n-i,n ^ 1- Recall that by Theorem 11.51 if y ^ N, then 
|G2n-i| — *■ 00 as n — *• cxo. 

Assume that p = y G N. In this case by Theorem 11.51 G'2n+i = lp,2n+i = iff g G 
Lg^g„(0, 1). Thus the statement i) has been proved. 

(ii) The function f^ satisfies the Eq. — /" + qf^ = Enfu on the interval [0, 1] and 
/a;(0) = /w(l) = ^,f'^^{lY = e^^" > for some /i„ G R. Each constant /i„ is so-called 
norming constant for the Sturm-Liouville problem —y" + qy = Xy on the interval [0, 1] with 
the Dirichlet boundary conditions y{0) = y{l) = [PT]. 

Recall that /U„,n ^ 1, is the Dirichlet spectrum of the problem —y" + qy = \y,y{0) = 
y{l) = on the unit interval [0, 1] and /i„ = En, n ^ 1. 



22 



Recall the well know result from [PT]: the mapping ^ '■ Q ^ [qo, (^n(Q'))i°5 (^«(?))i° ) is 

a real-analytic isomorphism between L^(0, 1) and R x IC x if. This gives the statement ii). 
(iii) If potential q is even, i.e., q G Lg„g„(0, 1), then each /i„ = 0, n ^ 1 (see [PT]). 

Recall the well know result from [PT]: the mapping ^e '■ q ^>- \ qo, {^n{q))i' ) is a real- 
analytic isomorphism between Lgj,g„(0, 1) and M x /C. This gives the statement iii). ■ 
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